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CALCULATION  OF  THE  AERODYNAMIC 
CHARACTERISTICS  OF  A  TWISTED 
JET 


R.  B.  Akhmedov  and  T.  B.  Balagula 

An  analysis  of  the  experimental  works  on  a  twisted  jet  and 
an  observation  of  various  theoretical  schemes  have  shown  that 
the  known  calculation  methods  permit  one  to  describe  a  twisted 
Jet  in  the  areas  where  the  initial  geometry  of  a  swirl  vane 
has  no  effect.  Meanwhile  the  area  of  a  flow  of  a  twisted  jet 
near  the  nozzle  depends  considerably  on  the  initial  conditions. 
Based  on  the  analytical  calculation  methods,  which  permit  the 
consideration  of  the  geometry  of  a  swirl  vane,  of  primary  interest 
is  the  method  of  equivalent  problem  of  the  theory  of  heat 
conductivity,  developed  by  L.  A.  Vulis.  The  applicability  of 
this  method  was  confirmed  for  a  whole  number  of  direct-stream 
fluxes;  however,  the  data  which  could  be  used  to  evaluate  its 
applicability  in  the  calculation  of  twisted  Jets  were  plainly 
Insufficient.  The  method  of  equivalent  problem  of  the  theory 
of  heat  conductivity  was  first  used  in  work  [5]  for  calculating 
velocity  fields  of  a  single  twisted  jet  flowing  Irom  a  helical 
tangential  twist  vane  with  the  twist  parameter  value  of  2.07 
In  this  case,  instead  of  the  usually  used  value  of  dynamic 
pressure  pu  ,  the  value  of  a  total,  pulse  was  introduced  into 
the  calculation  scheme  of  the  method,  which  included  the  excess 
static  pressure  (H«p„T+pu  ),  determined  by  the  equation  of  heat 
conductivity 
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(1) 


which,  with  an  arbitrary  profile  H  given  in  th$  form  H«Hn(y)  and 

3H  3  u 

boundary  condition  ^«0  with  y*0  and  y+O  with  y*«°,  permits  a 

solution  in  the  form 


The  basis  for  the  introduction  of  H  into  the  calculation 
scheme  of  the  method  of  the  equivalent  problem  is  in  the  fact 
that,  in  the  case  of  a  nonisobarlc  flow  (and  it  is  precisely  to 
this  class  of  fluxes  that  the  twisted  jets  belong),  it  is  not 

-j  ? 

the  flow  of  pulse  pu1*  that  is  the  remaining  value,  but  \  Hr  <lr — 

■  const,  i.e.,  H  is  the  remaining  value  in  all  cross  sections 
of  the  Jet. 

The  calculation  results  for  field  H  according  to  (2)  and 
H  values  obtained  by  the  test  data  [3]  have  shown  a  satisfactory 
agreement  for  a  twisted  Jet  flowing  from  a  helical  tangential 
swirl  vane  with  a  single  value  of  the  twist  parameter.  For  a 
more  extensive  verification  of  the  applicability  of  this  method, 
a  study  was  made  of  the  possibility  of  calculating  these  Jets 
using  the  various  techniques  for  twisting  these  jets. 

~  This  work  presents  the  results  of  the  experimental  and 
calculation  study  of  a  twisted  Jet  flowing  from  an  axial-flow 
swirl  vane  of  the  blade  type  and  from  an  axial-tangential  blade 
and  helical  swirl  vanes.  We  present  the  basic  structural  param¬ 
eters  of  the  studied  swirl  vanes  and  the  value  of  twist  parameters 
n,  determined  by  the  formulas  of  work  [1]. 


1)  Axially  bladed  register  A  (diameter,  in  mm,  of  neck  (d) 
200,  bushing  (dQ)--50,  dQ/d-0.25,  blades-16) 

blade  setting  angle  a,  in  degrees  Intensity  of  twist  n 


20 

30 

40 

50 

60 


0.35 

0.55 

0.79 

1.13 

1.65 


2)  Axially  tangential  blade  register  AT  (diameter  mm:  neck 
(d)-200,  bushing  (dQ)-70j  dQ/d~0.35) 

Angle  between  the  outlet 
blade  edge  and  the  axis 
of  the  cylindrical  channel 

3,  in  degrees  25  25  25  45  45  45 

Angle  between  the  blade 
and  the  tangent  to  the 


Internal  periphery  of 

twisting,  a,  in  degrees 

45 

30 

20 

45 

30 

20 

Intensity  of  twisting  n 

0.38 

0.64  0. 

80  0. 

63 

0.93 

1.35 

3)  Helical  register 

Y  (neck 

diameter 

d-200 

mm) 

iA/idth  of  tangential 

inlet  a,  in  mm 

128 

85 

128 

66 

100 

length  of  tangential 
inlet  b,  in  mm 

512 

340 

256 

265 

200 

Intensity  of  twisting  n 

1 

2 

2 

3 

3 

a/a 

0.64 

0.425 

0.64 

0.33 

0.50 

b/d 

2.56 

1.70 

1.28 

1.325 

1.0 

a/b 

0.25 

0.25 

0.50 

0.25 

0.50 
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The  calculation  for  the  fields  of  total  pulse  H  was  made. 

Just  as  in  [5],  according  to  formula  (2)  based  on  the  experimentally 
obtained  initial  profiles  of  complete  pulse  H.  In  accordance  with 
the  general  method  of  the  equivalent  problem  of  the  theory  of 
heat  conductivity,  relationship  5(x)  is  determined  by  comparing 
the  theoretical  and  experimental  data  with  a  certain  value  of 
the  transverse  coordinate  y.  If  this  comparison  is  usually  made 
when  y»0,  i.e.,  on  the  Jet  axis  then,  in  this  work,  relationship 
CH(x)  is  determined  according  to  the  coincidence  of  the  values  of 
positive  maximums  of  the  experimental  and  theoretical  profiles 
of  H  (Figure  1).  It  turned  out  that,  assuming  the  Invariance  of 
the  second  coordinate  n*y,  the  calculated  curves  have  shown 
regular  shift  in  the  calculation  maximum  ffom  the  experimental 
toward  the  area  away  from  the  Jet  axis  on  the  sections  close  to 
the  nozzle  (Figure  2)  in  the  case  of  strong  twisting  of  the  jet. 

The  following  transformation  was  introduced  to  calculate  this 
regular  displacement: 


i|  u  I-  a. 


(3) 


Value  a  is  selected  such  as  to  combine  the  calculation  and  this 
experimental  maximums . 


Figure  1.  Dependence 
✓Sjj-f(x)  for  a  twisted 

Jet  flowing  from  blade 
register  A.  Value  of 
angle  a,  in  degrees: 

1  -  20;  2  -  30;  3  -  40; 
4  -  50;  5  -  60. 


Figure  3  shows  the 
form  of  fields  H  in  an 
axial  register  after 
transformation  (3); 
experimental  fields  H 
have  two  special  features 


Figure  2.  Field  of  total  pulse  H  in  the  twisted  jet  flowing  from 
register  A  at  cx^^O0  .  1  —  experimental  data;  2  —  calculated,  using 

the  method  of  an  equivalent  problem  without  introducing  the 
transformation  of  the  displacement. 


■W  -05  II  Oi  1,0  to  3J0  -I  0  S3  Iff  1,5  -0.5  0 


Figure  3.  Field  of  total  pulse  H  in  a  twisted  jet  flowing  from 
register  A.  Values  of  a,  in  degrees:  a)  40,  b)  30,  1  -  test, 

2  -  theoretical  data. 


first,  maximum  H  is  displaced  away  from  the  axis  as  it  moves  away 
from  the  nozzle  and,  with  an  increase  in  twist  parameter  n,  this 
displacement  increases;  second,  with  large  values  of  the  twist 
parameter  the  H  retains  a  certain  constant  negative  value  in  the 
near-axis  region  in  each  section,  which  diminishes  in  absolute 


wwm».''y.w  -y . 


p  ../yffn  ^y< '  -VT  -T I TJ*  r>-  S  ?*  VK\y9® 


values  In  proportion  to  the  distance  from  the  opening, 


From  Figure  3  it  is  evident  that  the  calculation  curves  in 
the  area  of  positive  values  of  H  are  in  relatively  good  agreement 
with  the  experimental  data.  In  the  area  of  the  negative  H  values, 
the  agreement  between  the  theory  and  experiment  was  not  obtained 
even  after  the  introduction  of  transformation  (3).  Figure  4  shows 
the  calculation  results  for  the  H  fields  for  helical  and  axial-  ’ 
tangential  swirl  vanes  [2],  Just  as  in  the  case  of  art  axial  feed 
the  coincidence  in  the  near-axial  zone  cannot  be  considered 
satisfactory . 


An  analysis  of  the  obtained  materials  shows  that  the  nature 
of  rotation  of  the  air  flow  has  a  considerable  effect  on  the 
coincidence  of  the  theoretical  and  experimental  data:  a  better 
coincidence  is  obtained  for  a  helical  register  whose  field  of 
total  pulse,  especially  at  small  twistings,  is  close  in  form 
to  field  H  for  the  axial-flow  blade  register.  It  is  interesting 
to  note  that  in  the  volute  maximum  H  in  the  experiment  is  displaced 
to  the  jet  axis  (despite  the  fact  that  the  twist  is  relatively 
large,  n=l) .  In  this  case  we  did  not  have  to  Introduce  transforma¬ 
tion  (3).  Evidently  this  nature  of  field  K  in  the  helical  register 
can  explain  the  satisfactory  agreement  between  the  theory  and 
the  experiment,  obtained  in  work  (.5],  When  n=l  the  coincidence 
between  the  theory  and  experiment  in  the  region  close  to  the 
axis  is  worse  (Figure  4,  1,  a).  Such  a  divergence  of  data  with 
work  [5]  can  be  explained  by  the  fact  that,  in  addition  to  n  the 
nature  of  the  flow  field  (quasi-solid  or  potential)  is  affected 
oy  tne  value  of  the  ratio  of  tne  wiatn  of  tne  inlet 
connection  across  the  axis,  of  the  burner  to  the  length  of  the 
inlet  connection  along -the  axis  of,  the  burner  (|).  The  region  of 
quasi-3olid  rotation  which  is  characteristic  for  the  axial-type 
register  with  the  blade  pitch  angles  a>40°  increases  with  a 
decrease  in  g,  when  one  cannot  compare  the  results  of  the 
calculation  and  experiment  in  the  near-axial  region  also.  Since 
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Figure  4.  Field  of  total  pulse  H  In  a  twisted  jet  flowing:  ] 

from  a  helical  swirl  vane  [|-0.25;  n-1  (a),  n-2  (b)],  II  -  fr< 

the  axially  tangential  blade  swirl  vane  [a«30° ;  0-25° \  n-0 . 6 
(a);  a*45°,  n-0. 93  (b)3;  1  -  test,  2  -  theorectical  data. 

a  work  [5]  §»0.90,  while  In  our  work  g-0.25,  the  divergence  in 
he  results  can  be  understood  on  the  basis  of  these  concepts- 

An  analysis  of  the  data  shown  in  Figures  2-4  showed  that 
he  calculation  of  the  fields  of  total  pulse  H,  according  to  the 
ethod  of  unequivalent  problem,  yields  a  satisfactory  agreement 
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with  a  constant  transverse  coordinate  n*y  on  small  twists  of 

the  air  flow  (a*20  and  30°  In  register  A,  g«*0.25,  and  n*l  in  register 

U).  However,  the  calculation  of  a  strongly  twisted  Jet  (beginning 

with  angles  a>40°  in  register  A  when  !«0.25,  and  n«2  in  register 

U:  at  all  angles  a  and  3  in  register  A  B)  requires  the  introduction 

of  coordinate  y  into  the  theory  of  transformation.  The  inclusion 

of  static  pressure  pQT  into  the  calculation  scheme  of  the  method 

snould  have  taken  into  account  the  inverse  flow  of  the  fluid 

to  the  Jet  nucleus,  since  both  the  flow  and  rarefaction  on  the 

Jet  axis  develop  due  to  the  centrifugal  effect  in  the  twisted 

Jet.  However,  in  the  central  portion  of  the  Jet  we  are  unable 

to  obtain  an  agreement  between  the  calculation  and  test  data; 

consequently,  the  Introduction  of  value  H  instead  of  the  dynamic 
2 

pressure  pu  does  not  fully  take  into  account  strong  Inverse 
fluxes  arising  in  a  strongly  twisted  Jet. 

* 

The  calculation  of  the  field  of  axial  velocities  in  a 
twisted  Jet,  according  to  the  method  of  the  equivalent  problem 
of  the  theory  of  thermal  conductivity,  cannot  be  carried  out 
without  calculating  the  tangential  velocity  component,  since  the 


value  of  residual  static  pressure,  included  in  H»pcT+pu‘ 
be  determined  by  integrating  the  equation 


must 


ffV r_ 


(4) 


according  to  the  determined  relationship  u^(x,  y). 


As  is  known  the  tangential  velocity  in  a  twisted  Jet 
attenuates  considerably  faster  than  the  axial  velocity  (for  example, 
in  the  Jet-source  the  axial  velocity  decreases  with  distance  as 
j-  and  the  tangential  velocity  -  as  j2).  Therefore,  the  calcula¬ 
tion  methods  which  permit  one  to  obtain  a  continuous  deformation 
of  the  initial  profile  at  small  distances  from  the  opening 
acquire  a  special  Importance  in  calculating  the  tangential 
velocity . 


XlN 


Work  [53  purposes  the  use  of  the  equation  of  thermal  conduc- 

2 

tivity  for  value  u^  for  calculating  tangential  velocity  pu^.  The 
equation  analogous  to  (1)  is  written  for  this  equation 

«'*({>«•)  |  </((•«•) 

'  <*5*  ^  v  *  ,  (5) 

which  is  solved  with  the  following  initial  and  boundary  conditions: 


6t>«) 


it  «  /(?/) 

v  <>  i»«!  o. 

»-*•  °°  —y - *>0- 


2 

Boundary  condition  pu^«0,  when  y«0,  is  found  to  be  in  contradiction 

with  the  initial  equation  (5),  although  it  is  physically  valid. 

Actually,  when  y-0  the  solution  for  puf,  similar’ to  (2)  transforms 

9 


into 


•e  *  ^ 

1  •» 

and  is  never  equal  to  0. 


(7) 


Since  with  such  a  calculation  system  it  is  difficult  to 
assign  a  boundary  condition  where  u^-Q  on  the  jet  axis,  in  this 
work  we  purpose  a  new  method  for  calculating  the  tangential 
velocity.  The  hypothesis  lying  at  the  base  of  this  method  consists 


of  the  following:  the  equation  of  diffusion  is  satisfied  not 


by  value  pu 


2 

9* 


but  by  x  which  is  the  vector  component  of  vortex 


,  determined  by  the  formula 


"WTf‘iyvini!A'.w!f»"«’WT  'w|iii|,'’lj?a.«ff'1'  »■  ""TT?'  V.W  ''TPS?# 


In  the  hydrodynamics  of  a  viscous  laminar  motion  It. Is 
emphatically  proven  [4]  that  vortex  *«curl  v  satisfies  the  non- 
stationary  equation  of  thermal  conductivity,  which  for  the  x- 
oomponent  of  vortex  |  Is  written  In  the  form: 


vA*n 


where  v  -  viscosity;  A  -  Laplace  operator. 

Proceeding  from  nonstationary  equation  of  diffusion  (9)  to 
the  equation  of  the  method  of  the  equivalent  problem  using  the 
system  usual  for  this  method,  we  write  the  following  equation 
for  the  x-component  of  the  vortex 


<iif  i  i  rfiy 

-3r~”^r+7,“sr 


(10) 


The  initial  and  boundary  conditions  for  (10)  are  given  in  the  form: 


(*-0)  s  =»*,(?); 

S,>0  (*>0)  *  “►  0  when  y  -*•  oo. 


(ID 


Prom  the  connection  of  z  with  u.,  according  to  (8),  it  follows 

that  now  we  can  impose  the  necessary  additional 

condition  u.  on  u./„ _n»0.  Actually,  from  (8) 

9  9  y»u 


U,-Tl,yiv' 


(12) 


When  y-*-0 


-«o. 


Then  the  solution  of  equation  (10)  can  be  written  in  the  form 
similar  to  (2): 
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The  calculation  of  according  to  (13)  la  carried  out  numerically 
by  approximating  not  the  entire  Integrand  function,  but  only 
the  initial  profile  zQ(r).  This  permitted  us  to  reduce  the 
problem  of  calculation  of  to  calculating  with  the  aid  of  the 
so-called  P-functions  (5)*  The  calculation  data  and  their  comparison 
with  the  experiment  for  a  Jet  flowing  from  the  axial  blade 
apparatus  are  shown  in  Figure  5*  Profiles  z  were  determined  by 
the  formula 
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based  on  the  experimentally  measured  profiles  of  tangential 
velocity  u^.  It  is  evident  from  the  figure  that  based  on  the 
proposed  hypothesis  the  experimental  data  are  described  satis¬ 
factorily  by  the  theoretical  dependence. 

The  determination  of  vortex  z  by  the  experimental  profiles 
u^  in  the  opening,  in  a  number  of  cases,  leads  to  considerable 
errors.  The  second  term  (i/~)  included  in  expression  ( 1  h ) 
introduces  errors  into  the  initial  profile  z(0),  since  angular 
velocity  u  near  the  nozzle's  edge  drops  sharply  to  0  on  a  short 
segment  of  the  transverse  coordinate  y.  In  thooe  cases  when 
w  near  the  nozzle’s  edge  changes  slowly  (small  3lope  angles  of 
blades  in  a  twisted  Jet)  (Figure  6,a)  the  errors  in  profile 
z(0)  are  small.  With  large  slope  angles  (a>40-45° ) ,  w  changes 
sharply  near  the  nozzle’s  edge  and  the  coincidence  between  the 
calculation  and  test  worsens.  To  simplify  the  calculation  of 
u^  we  attempted  to  apply  a  new  method,  not  to  the  value  of 
vortex  z  according  to  (14),  but  only  to  the  value  of  angular 


,.J '•’l"  “M'p 


MCT^a.iaiiffTOia'awi 


Figure  5.  Field  of  the  x-component  of  vortex 
z  in  a  twisted  Jet.  Values  u,  in  degrees: 
a)  20,  b)  40;  1  -  test,  2  -  theoretical  data. 

UA 

velocity  “>■  y*  All  transverse  profiles  of  tangential  velocity  were 

recalculated  into  the  profiles  of  angular  velocity  according  to 

this  relationship.  The  initial  profile  of  the  latter  was  used  in 

(2)  to  calculate  the  angular  velocity  in  other  sections.  The 

calculation  results  are  in  good  agreement  with  the  experiment  (see 

Figure  6).  This  method  of  calculation  is  not  sufficiently  strict, 

00 

since  u  is  not  a  constant  value;  however  j  w  ydy  maintains  a  constant 

(with  an  accuracy  up  to  10-2U£)  value  in  all  these  sections.  In 
practice,  to  find  u^,  using  the  angular  velocity,  is  simpler 
and  this  method  can  be  recommended  for  engineering  calculations. 

Dependence  ^^(a),  necessary  for  the  calculation  is  shown  in 
Figure  7.  From  the  fulfilment  of  equation  (4)  (the  results  of 
its  experimental  verification  are  shown  in  Figure  8)  stems  the 
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Figure  6.  Field  of  angular  velocity 
u  at  the  outlet  from  register  A. 
Values  a,  in  degrees:  a)  40,  b)  20; 
1  -  test,  2  -  theoretical  data. 


possibility  for  calculating  the  field  of  static  pressure  p  in 
a  twisted  Jet. 


Figure  7.  Dependence  /T"(x)  for  a  twisted 

Jet  from  blade  register  A.  Values  a,  in 
degrees:.  1  -  20,  2  -  30,  3  -  40. 


After  calculating  u  according  to 
the  proposed  method  and  field  H  according 
to  (2),  it  is  possible  to  calculate  the 
fields  of  the  axial  velocity  component 
and  to  complete,  in  this  manner,  the 
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calculation  of  the  principal  aerodynamic  characteristics  of  a 
twisted  Jet. 


Pig.  8.  Experimental  varification  of 
the  equation  (4)  for  a  twisted  Jet 
flowing  from  register  A  with  a*»40 
(section  x-1.0) . 
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